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Abstract 

In this paper, the synchronization problem for a class of discrete-time complex-valued neural networks with time-varying 
delays is investigated. Compared with the previous work, the time delay and parameters are assumed to be time-varying. By 
separating the real part and imaginary part, the discrete-time model of complex-valued neural networks is derived. 
Moreover, by using the complex-valued Lyapunov-Krasovskii functional method and linear matrix inequality as tools, 
sufficient conditions of the synchronization stability are obtained. In numerical simulation, examples are presented to show 
the effectiveness of our method. 
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Introduction 

In the past decades, the study of dynamical neural network has 
attracted researchers' attentions because of its potential applica- 
tions in a variety of areas, such as image processing, combinatorial 
optimization problems, pattern recognition, signal processing and 
so on (see, for instance [1-5]). In the investigations of dynamical 
neural network, synchronization, control and stability analysis 
have attracted people's interests in the research of chaotic systems, 
complex nonlinear systems and dynamical neural networks (see 
[6-1 1] and references therein) since Pecora and Carroll achieved 
synchronization between two chaotic oscillators by PC (Pecora 
and Carroll) method [12]. As we know, time delays commonly 
exist in the neural networks because of the network traffic 
congestions and the finite speed of information transmission in 
networks. So the study of dynamic properties with time delay is of 
great significance and importance. However, most of the studied 
networks are real number valued. Currendy, in order to 
investigate the complex properties in complex-valued neural 
networks, some complex-valued network models are proposed. 
For example, Hu gives the global stability of complex-valued 
recurrent continuous neural networks with time-delays [13] and 
Zhou further studies the boundedness and complete stability of 
complex-valued neural networks [14] in 2013. For discrete-time 
complex-valued neural networks, the boundness and stability of 
neural networks without time delays [15,16] and the boundness 
and stability of neural networks with time delays are studied by 
Zhou and Duan. However, in aforementioned discrete-time 
complex-valued neural networks, time delays are fixed and 
parameters are constant, which are not suitable. It should be 



noticed that time delay is a common phenomenon in network 
because of the signal transmission. Generally, time delays can be 
divided into time-varying time delays and constant time delays. 
Time-varying time delays are difFerent from constant delays 
because the delay varies with time. In practice, time delays always 
vary in a bounded range instead of a fixed point and the constant 
delays can be seen as a special case of time-varying delays. So the 
study of time-varying delays has more potential applications. On 
the other hand, the study of time-varying delays is more 
challenging, the derived theorems of time-varying delays can be 
easily applied in the network with constant delays when upper 
limit equals to the lower limit as T = T. However, the corresponding 

theorems of constant delays are hard to be used in the study of 
time-varying delay network. 

Motived by above discussions, by separating the real part and 
imaginary part of the neural networks and constructing complex 
Lyapunov-Krasovskii functional candidates, we will investigate the 
synchronization problem for a class of discrete-time complex-valued 
neural networks with time-varying delays. Although some sufficient 
conditions for stability and synchronization problems of discrete real 
number neural networks have been derived by some researchers, as 
far as we know, there has been no literatures investigate the discrete- 
time complex-valued neural networks with time-varying delays. We 
believe that the synchronization problems for such kind of networks 
are still remains open and challenging. 

The rest paper is organized as follows. In section 2, the basic models, 
preliminaries and lemmas are presented. Section 3 presents stability 
analysis and sufficient conditions with linear matrix inequality (LMT). 
Some numerical simulations and examples to show the robustness and 
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effectiveness of our methods are in section 4. Finally, we give some 
concluding remarks in section 5. 

Notations: Throughout the paper, R" represents the n-dimension- 
al Euclidean space. R nxm is the set ofnxm real matrices. T means 
the transpose of the corresponding matrix and the symmetric 
matrix X>0 (respectively, X<0,X>0 and X<0) means that X 
is positive semidefinite (respectively, negative semidefinite, positive 
definite and negative definite). diag{ - ■ •} denotes a block-diagonal 
matrix and • is used to represent a term induced by symmetry. If 
not explicitly stated, matrices dimensions are assumed to be 
compatible for algebraic operations. 

The system model and preliminaries 

In this paper, we will consider the following discrete-time 
complex-valued neural network model consisting of n coupled 
nodes with time-varying delays: 

x(k+ l) = A<j(x(k)) + B<j(x(k-x(k)))+H, (1) 

where x = [x\,xi, ■ ■ ■ ,x n ) T k the state vector, wis the number of 
neural cells. A=(a») eR" xn , B=(b ij ) eR nx " are the 

V V I nxn ' V 'J/nxn 

connection weight matrix and the delayed connection weight 

matrix. H= (h\,li2, ■ ■ ■ ,h n ) T represents the input vector, 
cr(x) = max{0,Re(x)}+j-max{0,Im(x)} is a complex-valued 
function, j = \J — 1. Time delay x(k) ranges from x to T as 

x <%{k)<%. Complex-valued parameters in the neural network 
can be represented as a i j = a R +ja\-, h„ b R ■ //>'. Then the 
model can be separated as 

xf{k+ 1)= £ 4°(x?(k))- Y, 

7=1 7=1 

+ Ybf ja (xf(k-x(k))) (l_l) 

7 = 1 

-Yb\M(k-<k))) + hf, 

7 = 1 

x\(k+\)=Y <$*tf(k)) + Y 4< x f( k )) 

7=1 7=1 

+ Ybfc(x' j (k-T(k))) (1-2) 

7=1 

+ Yb' ij o{xf{k-x{k)))+h I i , 

7=1 

where xf and x\ are the real part and imaginary part of variable 
Xi, respectively, afj and are the real part and imaginary part of 
connection weight ay and bfj and by are the real part and 
imaginary part of delayed connection weight by. hf and h\ arc the 
real part and imaginary part of input hj. Connection weight 

matrices are represented as A R =[af j \ eR" xn , 
^=(4) eR" xn , B R =(bf) eR nxn and 

\ • / nxn \ • / nxn 

B 1 = \b\j\ eR" x ". Consider two discrete-time complex-valued 

neural networks with time-varying delays labeled as (p and 
Variables at step k in these two neural networks are defined as 
x(k,(p) and x(k,\jj) and x R (k,(p), x R {k,\j/), x'(k,<p) and x 1 {k,\\i) are 



corresponding real parts and imaginary parts. It follows from the 
models (1-1) and (1-2) that 

x R (k+\,<p)-x R (k+\,\jj) 
= A R (o(x R (k, ( p))-cj{x R {k,m 

-A'ioix'ik^-aix'ikM) (2 - 1) 

+ B R (a(x R (k - x{k),<p)) - a(x R (k - x(k)M) 

- B'{a{x'{k - %(k),<p)) - a(x'(k - x{k),xj/))), 

x r (k+l,(p)-x J (k+l,<l/) 
= A R (<j(x'(k,(p)) - <j(x r (k^))) 

-A r (a(x R (k,<p))-<j(x R (kM) (2-2) 
+ B R (a(x'(k - x(k),<p)) - aix'ik - x(k),\ji))) 

- B'{a{x R {k - x(k),<p)) - a(x R (k - x(k),\jj))), 

where y R (k) = x R (k,(p)-x R (k^), y'(k) = x'(k,(p)-x'(k^), 
f R (k) = o(x R {k,<p)) - a(x R (kM, f\k) = G{x'(k,<p)) - o(x\kM- 
Then we have 

y R (k+l) = A R f R (k)-A I f(k) 

+ B R f R (k - x(k)) - B'f'{k - x(k)) 

(3) 

y'(k+l) = A R f'(k) + A I f R (k) 

+ B R f I {k-x(k)) + B I f\k-x{k)). 

The initial condition associated with complex-valued network 
(3) is defined as 

Xj(s) = (f>i(s),se[ — x,0] ,i = 1 ,2, • • • ,n, 

where Re(i^,(.y)) and Im(<^ ,(.?)) are continuous when xe[ — x,0]. 

In this paper, following definition and lemmas are listed to get 
the stability conditions. 

Definition [16]: 

The equilibrium point x of the model (1) with the initial 
condition Xi(s) = <l>i(s) is said to be globally exponentially stable if 
there exist two positive constants M. > 0 and 1 > £ > 0 satisfies 

\\x(k,(p) — x\\<Me k sup \\<p— x\\. 

.se[-t,0] 

Lemma 1 (Schur complement, [17]): 
Given the following matrix 

27= El ^ 

Zi -Zi _' 

where 27iis a non-singular matrix and S\ =S\, Z2 >0 and £3 is a 
constant matrix, then we say Ei+SlZ^^S^ is the schur 
complement of 27 about Z\ and we have the following conclusion: 
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holds if and only if the following schur complement holds 

27<0 



Lemma 2 [18]: 

Consider a symmetric positive-semidefmite matrix >PeR n><n 
(that is to say, f T = f > 0), scalar a, > 0(7 = 1 ,2, • • • ) and vector 
XiGR". We can get the following inequality 



+ O0 +O0 + GO +CO 

( ^ a iXi ) T W{ a,x,) < ( «0 E a <*7 y *<- 



Stability for Discrete Complex Neural Networks 

Proof: The detailed proofs for eq. (4) can be found in 
the Appendix SI. 

As we mentioned in the introduction, time-varying delays are 
completely different from constant delays. Constant delays can be 
seen a special case of time-varying delays when the delays vary in a 
fixed point. On the other hand, some existing methods can not be 
applied in time-varying systems. In order to show the difference, 
the following corollary will be presented. 

Corollary 1 

For a constant time delay T, the dynamical neural networks (1) 
will be globally stable if there exist positive matrices P\ > 0, 
P 2 >0, Qi>0,Q 2 >0, H\ >0 and TZ 2 >0, diagonal matrices 
Pi >0, © 2 >0, £i >0, and C 2 >0, such that the following LMI 
holds. 



i = i ;=i 



Assumption [13]: 

gi(»),i= 1,2, • • • ,n satisfy the Lipschitz continuity condition in 
the complex domain, which means there exists a positive constant 
E/,i= 1,2, • • • ,n for any state variable u and v, one has 



\gi(u)-gi(v)\ <ei\u- v\,i= 1,2, •• • ,«. 



Where £; is the corresponding Lipschitz constant. 

Methods 

In this section, we will deal with the synchronization problem of 
the aforementioned complex-valued neural networks (1). First, we 
will give the main result in this paper as follows. 

Theorem 1 

The dynamical neural networks (1) will be globally stable if 
there exist positive matrices P\ >0, 7 :> 2>0, Qi>0 and Q 2 >0, 
diagonal matrices V\ > 0, T>2 > 0, £\ > 0, and £2 > 0, such that the 
following LMI holds. 



-V + (l+t-T)Q + V 



0 

Q + £ 



Q 



0 
0 

V 0 
0 £ 



. (4) 



-v+Q+{t-\)n+v 



-K 
x 



-TZ 



0 



V 



(V 0\ 



(5) 



where *¥ - 



n 



-A R Ri 
-A'R 2 



-B R R 1 
-B'R 2 



B'Ri 

-B R R 2 



and 



Other variables and matrices are the same 



A'Ri 
-A R R 2 

TZi 0 

0 ft 2 

with definitions in theorem 1. 

Proof: For the stability of a complex-valued neural network 
with constant time delay T, with the similar method in [16], it is 
easy to get the corollary. 

Remark 1 

The class of complex-valued neural network with constant 
time delay has been investigated in [16] and theorem has been 
derived. Different from the previous study, the corollary 1 is 
established by separating the real part and imaginary part from 
the complex-valued neural network. Which means the problem 
of complex-valued neural network is translated into the 
stability of corresponding real neural networks and the 
network data are easier to handle. However, in practice, time 
delays are always time-varying, and then the above corollary 
will be infeasible. 

Above linear matrix inequalities are established under the 
circumstance of not taking the assumption into account. If above 
assumption, which is an extension of the real-valued function, is 
satisfied. We have 



where 



V = 



'V\ 0 

, 0 v 2 



0 



0 

p 2 



Ci 0 

0 £ 2 



Q = 



Qi 



0 

Q2 



and 



Corollary 2 

The dynamical neural networks (1) will be globally stable if 
there exist positive matrices V\>0, 7 5 2>0, Q\>0 and Q2>0, 
diagonal matrices 2?i>0, T>2>0, £i>0, £2>0, positive 
constants £ > 0, such that the following LMI holds. 



Q-- 



/ A R V\A R + A 1 V 2 A I -A R V\A' +A I V 2 A R A R P l B R + A I P 2 B I -A R V\B' +A I P 2 B R \ 

• A'PxA 1 +A R P 2 A R -A'P l B R + A R P 2 B I A'PiB 1 + A R P 2 B R 

• • B R P X B R + B'P 2 B' -B r P x B' +B'P 2 B r 

• • • B'PiB 1 + B R P 2 B R ) 
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4> = 



/ -P + (1+Z-t)Q + £ 2 V 



\ 



o 



n- 



o 

0 

V 0 

0 c 



,(6) 



where V, Q, V, C and Q have the same meaning as in theorem 1. 

Proof: The detailed proofs for eq. (6) can be found in 
the Appendix S2. 

Remark 2 

Compare the theorem 1 with the corollary 1, we can find that 
the theorem is the special case when the Lipschitz constant e = 1 . 
For real-valued neural network, the activation function is always 
chosen smooth and bounded. However, in complex-valued 
network, according to Liouville's theorem [19], the activation 
functions cannot be both bounded and analytic. Which means 
theorem 1 is less conservative and more general in practice. 

Results and Discussion 

In this section, examples are provided to demonstrate the 
robustness and effective of our method. 

Example 1 

Consider a two-neuron complex- valued network, where 



A = 



B- 



0.2+0.1/ 
-0.1+0.1/ 



-0.2-0.2/ 
L 0.1+0.3/ 



-0.2 + 0.2/ 
0.2 



-0.1 
-0.1+0.2/ 



H=[6+j 3-2/] T 



The time delays are bounded as T = 1 and T = 2, then from 

theorem 1, by using the LMI toolbox, we can solve the LMI (3) 
with the proposed matrices and parameters and the feasible results 
are shown as 



Qi = 



0.6107 0.045 
0.045 0.5565 



0.598 0.0465 
0.0465 0.5512 



£>i=diag(0.5388,0.5014), 



£i=diag(0.4213,0.3648), 



D 2 = diag(0.5294,0.4968), 



£ 2 = diag(0.4137,0.3642) 



Set the initial states as x\\ = 8 + 4j,xn = 2 — 3j, xi\ = 2+j and 
X22=—5+j. Errors are defined as e\=x\\— xj\ and 
e 2 =-" : i2 — Based on the theorem 1, the discrete-time coupled 
neural networks with mixed time delays will get the synchroniza- 
tion and it is proved by numerical simulations. Fig 1 shows the real 
part states of discrete-time complex-valued neural networks and 
Fig 2 shows the imaginary part states of discrete-time complex- 
valued neural networks. Fig 3 depicts the errors of discrete-time 
complex-valued neural networks. 

Example 2 

Consider a neural network with three nodes, where 



A = 



0.1+0.2; 0.2-0.1/ 0.1 
0.3 + 0.2; —0.1—0.1/ 0.2 + 0.1/ 
0.1 -0.1—0.1/ 0.2 + 0.1/ 



v 2 



2.1093 0.1353 
0.1353 1.9460 



2.0665 0.14 
0.14 1.9258 



B- 



0.2/ 0.1+0.1/ 0.1+0.1/ 
0.2+0.2/ 0.1+0.1/ 0.2 + 0.1/ 
0.1-0.1/ -0.2/ 0.1+0.3/ 



H=[3-j -2/ 2 + 2/]' 
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Figure 1. Real part of the two-neuron complex-valued neural networks. 

doi:10.1371/journal.pone.0093838.g001 



Time delays are bounded as t = 2 and T = 4, then from theorem 

1 or corollary 2, by using the LMI toolbox, we can solve the LMI 
with the proposed matrices and parameters and the feasible results 
are shown as 



22.1626 -6.0156 2.0340 
Qi= -6.0156 25.6346 -4.5197 
2.0340 -4.5197 26.9903 



Pi 



114.7976 -27.1947 9.1416 
-27.1947 139.5118 -20.6598 
9.1416 -20.6598 99.7014 



17.2688 -0.2016 -1.0305" 
Q 2 = -0.2016 16.0091 -0.0880 
-1.0305 -0.0880 27.7315 



v 2 -- 



85.0556 -1.1663 -5.6413 
-1.1663 75.5833 -0.5345 
-5.6413 -0.5345 117.2383 



D,=diag(43.3732,42.676,15.5231), 



£i=diag(19.1035,12.7346,25.1787), 



5 
4 

M 3 

£ 2 
M 

1 

j| -1 
-2, 



- 




-»-lm(x11) 




7 


-^Im(x12) 






— Im(x21) 






— -Im(x22) 













10 



15 

k 



20 



25 



30 



Figure 2. Imaginary part of the two-neuron complex-valued neural networks. 

doi:1 0.1 371 /journal.pone.0093838.g002 
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Figure 3. Errors of the two-neuron complex-valued neural networks. 

doi:10.1371/journal.pone.0093838.g003 



2? 2 = diag(31.0615,26.3729,28.9058), 



£ 2 =diag(l 6.421,15.595,25.704) 



Conclusion 

In this article, we have presented a new approach to deal with 
the stability problem of discrete-time complex-valued neural 
networks with time-varying delays. Compared with existing 
results, the derived theorems and corollaries are easier to be 
solved with the Matlab LMI toolbox because real part and the 
imaginary part of complex-valued variables are separated and 
derived results are real number style. What's more, we propose a 
feasible method to solve time-varying delays of discrete-time 



complex- valued neural networks. Constant time delays are 
considered as a special case of time-varying delays in our paper 
and also solvable. In addition, it should be noted that time-varying 
delays are bounded and lager time-varying ranges will lead the 
LMI unsolvable. Finally, simulations show the effective and 
robustness of our method. 

Supporting Information 

Appendix SI Proof of eq. (4). 

(DOC) 

Appendix S2 Proof of eq. (6). 

(DOC) 
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